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Dark Matter particles with inelastic interactions are ubiquitous in extensions of the Standard 
Model, yet remain challenging to fully probe with existing strategies. We propose a series of powerful 
searches at hadron and lepton colliders that are sensitive to inelastic dark matter dynamics. In 
representative models featuring either a massive dark photon or a magnetic dipole interaction, 
we find that the LHC and BaBar could offer strong sensitivity to the thermal-relic dark matter 
parameter space for dark matter masses between ~ 100 MeV-100 GeV and fractional mass-splittings 
above the percent level; future searches at Belle II with a dedicated monophoton trigger could also 
offer sensitivity to thermal-relic scenarios with masses below a few GeV. Thermal scenarios with 
either larger masses or splittings are largely ruled out; lower masses remain viable yet may be 
accessible with other search strategies. 


I. INTRODUCTION 

The observed cosmic abundance of dark matter (DM) 
mm is clear evidence of physics beyond the Standard 
Model (SM). While the non-gravitational dynamics of 
DM are not currently known, additional interactions with 
the SM are well-motivated and arise in many theories. In 
thermal DM scenarios - a compelling paradigm for DM 
physics - these possible interactions account for the ob¬ 
served cosmological abundance via DM annihilation into 
the SM. This framework motivates DM physics with both 
mass scales and interaction strengths potentially accessi¬ 
ble at current experiments. It is therefore imperative to 
robustly test the thermal DM paradigm with broad and 
complementary experimental approaches. 

In the most commonly studied weakly interacting mas¬ 
sive particle (WIMP) DM scenario, there is only one DM 
particle which interacts with the SM via a single type 
of interaction. In this case, obtaining the observed DM 
abundance through thermal freeze-out fixes a minimum 
coupling between SM and DM states, and a variety of 
experiments can be used to test the possibility of ther¬ 
mal DM. These constraints are strongest for DM with 
masses above a few GeV, and a combination of direct- 
detection, indirect-detection, collider, beam-dump, as- 
trophysical, and cosmological probes can dramatically 
narrow the window for thermal DM EHB]. 

However, the dark side of particle physics could exhibit 
a richer structure, especially given the complexity of the 
SM, and DM could even live in a dark sector (DS) with 
additional particles and forces [MU- This presents both 
new challenges and new opportunities: some probes of 
DM can be dramatically less sensitive in even the sim¬ 
plest DS scenarios, relaxing the constraints on thermal 
DM, while at the same time new prospects for the dis¬ 
covery of DM emerge. 

In this paper, we explore some of the striking signa¬ 
tures at colliders that can appear in a generic DS. A 
representative example of a DS consists of a dark mat¬ 
ter particle which is charged under a hidden gauge or 
global symmetry. The DM can have both a symmetry¬ 
preserving mass and, if the symmetry is spontaneously 


broken, also a symmetry-violating mass, which splits 
the mass eigenstates. In the limit that the symmetry¬ 
breaking mass is much smaller than the symmetry¬ 
preserving mass, the DM interactions are off-diagonal 
(between different mass eigenstates). This is a straight¬ 
forward realization of the inelastic DM (iDM) scenario 
proposed by Tucker-Smith and Weiner [12] . with pro¬ 
found implications for experimental probes of DM. In 
particular, the abundance of the heavier eigenstate can 
be large in the early universe, facilitating efficient co¬ 
annihilation of DM, whereas the heavier eigenstate is de¬ 
pleted today, suppressing indirect- and direct-detection 
signatures. The small DM halo velocities imply that 
DM has insufficient energy to up-scatter into the heavier 
state, and so interactions through the off-diagonal cou¬ 
pling are ineffective today. 

By contrast, the energies of colliders such as the LHG 
and H-factories are typically large enough to produce 
both the lighter and heavier DM mass eigenstates. In the 
iDM scenario, the dominant DM coupling to SM states is 
through the off-diagonal interaction, and so both eigen¬ 
states are produced simultaneously. When the heavier, 
“excited” dark state (denoted with an asterisk) decays to 
the lighter, “ground” dark state, some visible SM states 
are emitted. Thus, in addition to the standard DM miss¬ 
ing transverse energy (Xt) collider signature [T3H39], 
where the DM system recoils off of a jet, photon, vec¬ 
tor boson, or Higgs boson, iDM models typically feature 
the emission of associated soft SM states HU HO]. The 
characteristic iDM collider signature is the production of 
DM -I- DM* in association with a hard SM object X, fol¬ 
lowed by the subsequent decay of DM* —)■ DM -|- Y for 
some potentially different SM states Y. The production 
is summarized as 

pp A-kDM-kDM* 

^ X + DM + ^DM* ^ DM -I- =X + ^t + Y , 

and is depicted schematically in Fig. X is any state 
that can be used to trigger on the event and reconstruct 
throughout this study, we consider the case where X 
is a jet for hadron colliders, and X is a photon for lepton 
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colliders. V depends on the mode by which DM couples 
inelastically to the SM. As we elaborate in Sec. |llj the 
representative models we consider lead to two promising 
modes of DM* decay, namely 1* = 7 and V = . 

In this study, we propose a suite of collider searches 
for inelastic DM signatures. In particular, we focus on 
DM and DM* masses in the 100 MeV-tens of GeV range, 
and splittings of order ~ 1 — 10% of the DM mass, one 
of the blind spots of the current search program due to 
the suppression of indirect and direct detection signa¬ 
tures. For such light masses, when the DM and DM* 
states recoil against a comparatively hard jet or photon, 
the soft SM decay products of the excited state are typ¬ 
ically aligned with the missing momentum. We show 
that this feature allows for the effective suppression of 
the electroweak backgrounds for conventional mono jet 
and monophoton searches. Moreover, for thermal-relic 
DM-SM couplings and 0(10%) mass splittings between 
the ground and excited states, the decay of DM* can also 
occur on macroscopic distances, leading to displaced ver¬ 
tices and other non-prompt phenomena. This results in 
the possibility of a low-background search over much of 
the DM parameter space; indeed, the distinctive kine¬ 
matics of iDM production at colliders allows for sensi¬ 
tivity to the interactions responsible for the cosmological 
DM abundance. This is in contrast with traditional col¬ 
lider probes of many elastic DM models, where the large 
SM backgrounds strongly limit the sensitivity to thermal 
relic scenarios where the dominant DM-SM interaction 
is mediated by a new particle with mass at or below the 
weak scale. We illustrate the sensitivity to iDM in two 
concrete representative models: a model where DM in¬ 
teracts with the SM via a kinetically mixed dark photon 
M: and a model where DM couples inelastically to the 
SM via a magnetic dipole moment |1H US] . We summa¬ 
rize our results in Figs. 

In addition to the iDM signatures considered here, 
there are many other manifestations of dark sector states 
at colliders, which can give taggable objects such as hard 
final-state radiation of new gauge bosons, energetic SM 
states from excited DM decay, and dark showers [iiHSD] 
that are complementary to our studies. Inelastic DM 
decays with monojet -I- soft hadronic displaced signa¬ 
tures were considered in the contact interaction limit in 
Ref. IH]. Finally, monojet -|- soft object searches are also 
useful for compressed supersymmetric spectra (for recent 
examples, see Ref. [BTHHO] !. although our work examines 
parametrically different masses and splittings and differ¬ 
ent final states, focusing particularly on long-lived decays 
and exploiting different kinematic features. 

The rest of this paper is organized as follows: in Sec.[TTl 
we present the two classes of representative models that 
this paper studies. We then propose a series of poten¬ 
tially powerful collider searches at both R-factories and 
the LHC in Se c.|III| The cosmology of these models is de¬ 
scribed in Sec. |IV| Finally, we discuss existing constraints 
on the simplified models in Sec. 



FIG. 1: Schematic diagram depicting a characteristic iDM 
production event at the LHC: pp ^ j + DM DM*, where DM 
is the dominant DM component in our halo and DM* is a 
heavier, unstable DS state. The final state contains visible SM 
particles (7 or l^£~) and missing transverse energy produced 
in association with a QCD jet. At lepton colliders, a similar 
process of interest is e^e“ —>-7 DM DM*. 


II. REPRESENTATIVE MODELS 

The classes of models we consider in this paper all fea¬ 
ture dark matter currents that couple inelastically to the 
SM. Some of the simplest realizations of DM-SM interac¬ 
tions involve an additional massive mediator particle that 
connects SM and DM currents - e.g. a kinetically mixed 
dark photon or a We also consider the possibility that 
DM couples to SM gauge bosons via higher-dimensional 
operators. In the simplest such example, fermionic DM 
couples to 7 /Z via a magnetic dipole moment: if DM is 
Majorana, this coupling must be inelastic. 

We now discuss in turn each of the two models in 
our study: a dark photon and magnetic inelastic DM 
(MiDM). 


A. Dark Photon Model 

A simple, well-motivated candidate mediator between 
the dark and visible sector is a massive dark-sector gauge 
boson A', whose most general renormalizable Lagrangian 
contains 

+ gB^J^ + goA'^J^. ( 1 ) 

Here, ey is the kinetic mixing parameter. A' is the mas¬ 
sive “dark photon” of a broken U{1)d symmetry, B is 
the hypercharge gauge boson, and B^^ are the dark- 
photon and hypercharge field strength tensors, Jy is the 
SM hypercharge current, Jb is the dark current, and raA' 
is the dark photon’s mass. 

After electroweak symmetry breaking, A mixes with 
both 7 and Z, so in the SM mass eigenbasis the La- 
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FIG. 2: Collider projections for fermionic iDM in the dark photon model with ao =0.1 and mA' Imi = 3 vs. thermal relic 
density target and other constraints. For LHC projections (red dashed), we consider a jet + + displaced lepton-jet topology 

in 13 TeV running with 300 fb“^. For B-factory projections, we consider existing constraints from BaBar on photon -|- (green 
solid), projected reach of photon -I- displace d lep ton signatures (green dashed), and projections for a possible Belle II 

we also show the projection for a proposed 
since this search would veto visible energy 


monophoton -|- search (purple dashed). See Sec. Ill for details. For A = O.lmi, 
fixed-target missing-momentum experiment (orange dashed) drawn from Ref. |61| 
from X 2 de-excitation, we conservatively assume it only has sensitivity to A = O.lmi. Also shown are constraints from LEP 
[H2] and {g - 2)^ [^, whose sensitivities do not scale with y\ see Sec.[^ Both experimental constraints are only sensitive to the 
visible coupling t and m^'. To avoid overstating these bounds, we conservatively show their y contours for the reasonably large 
values oi 0.0 and niA' lrn\ given above, which reveals most of the allowed parameter space (see Sec.|TI|. For smaller values of 
OLoijnx!raA')'^■, as shown in Fig. the y-reach for these bounds is greater and shifts linearly downwards to cover more of the 
thermal relic line. The jagged spikes represent annihilation to hadronic final states as discussed in Appendix 


grangian contains 

^ . ( 2 ) 

After diagonalizing the kinetic terms, the dark photon’s 
couplings to SM fermions are approximately given by [5B] 

m| cos OweQf - m\,gYf 


where {Yf)Qf is the SM fermion’s (hyper)charge. In the 
limit of a light A', the mixing is predominantly with the 
photon and QA'ff ^ ey cos ^vveQ/, so the visible sector 
acquires a millicharge under U{1)d and we exchange ey 
for the related parameter 

e=eycos6>w. (4) 

After U{\)]A symmetry breaking, the DM charge eigen¬ 
states will generically mix, giving rise to a split spectrum 
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Fermion Thermal Relic iDM, A = 0.1 nii, ao = qtem Fermion Thermal Relic iDM, A = 0.2 mi , ao= (Tem 





FIG. 3: Collider projections for fermionic iDM in the dark photon model with varying q_d and mA'/mi. The plots in this 
figure depict the same mass range and mass splittings as the top row of Fig. Top row: same /mi ratio as Fig. but 
with q;_d = o instead of an = 0.1. Bottom row: same otn as in Fig. but with the DM-mediator mass ratio niA'/mi = 10. 
Similar scaling applies to the scalar scenario shown in Fig. 


and inelastic DM, and we show the spectrum and inter¬ 
actions below. 

In addition to proposing collider searches for DM cou¬ 
pled via A', we also explore how the collider constraints 
compare to the parameters giving the observed relic 
abundance and other constraints. Since the A'-mediated 
scenario depends on five parameters - the lightest DM 
mass, TOi; the DM mass splitting. A; the A' mass, mA'] 
the dark gauge coupling, an = ffi)/47r; and e - care must 
be taken to avoid overstating bounds on the parameter 
space. In the A ^ mi < niA '/2 limit, the DM annihi¬ 
lation rate largely depends on only two parameters: mi, 
and the dimensionless interaction strength, y: 

2 f mi Y _ (.1 

av oce ao [ - ] =y , (5) 

\mA' J 


which is insensitive to individual choices for each param¬ 
eter so long as their product remains fixed [ST]. Small 
values of an or of mi/iriA' would lead to an overabun¬ 
dance of dark matter unless e is correspondingly larger; 
on the other hand, different experimental bounds may 
not scale straightforwardly with y. For example, preci¬ 
sion QED constraints depend only on e and m^', and are 
independent of mi and ao- These constraints, expressed 
in terms of y, would therefore be overstated for small val¬ 
ues of an and mi/niA' relative to the y value required 
for the observed relic abundance. To be conservative, it 
suffices to choose large, order-one values of these quan¬ 
tities in computing experimental bounds on y. We show 
later how the results scale for different values of an and 
mi/mA'- 

For the secluded DM scenario (m^' < mi) |5], the 
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FIG. 4: Same as Fig. but for scalar iDM in the dark photon model. 


annihilation process DM + DM A' + A' sets the relic 
abundance, which is independent of the A' coupling 
to SM states. Thus, there is no robust experimental 
target for this scenario^. It’s possible, however, to still 
produce DM + DM* through a virtual A' at colliders, 
with subsequent decay DM* —DM + A'. The A' 
subsequently decays into SM final states. This kind 
of topology would fall under the scenarios studied by 
Refs. [48l|49]. Furthermore, the direct production of the 
DM states may not be the discovery channel of this class 
of models, as now the A' could be produced directly and 
observed through its decays into or into dijets (see 
Refs. [Ml ES] for recent studies). Because ruA' ^ 2 toi 
offers a clear, experimentally promising target for the 


^ For a discussion of the interpolating regime mi < m^/ < 2mi, 
see Refs. 16711681 . 


parameters giving the observed relic abundance, we 
focus on that scenario. 


Returning to iDM, the y necessary for freeze-out grows 
with increasing mass splittings. A, but is still a useful 
variable to characterize the parameter space for fixed A. 
For a purely inelastic coupling, the A > mi regime is ex¬ 
cluded by a combination of collider and precision-QED 
probes (see Figs. pJQ. Similarly, for sufficiently small 
DM parameters m^JmA' or ^ 1, these same con¬ 
straints rule out the thermal freeze-out hypothesis (see 
Fig. §. Thus, the viable parameter space for thermal 
iDM coupled to an A! requires A < 0(10%), comparable 
DM/mediator masses, and sizeable «£> < 1, so our search 
strategy in this paper primarily targets this regime. 
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MiDM Thermal Relic, A = 0.03 mj MiDM Thermal Relic, A = 0.05 m\ 



FIG. 5: Collider projections for dipole IDM vs. thermal relic density target and other constraints. For LHC projections, we 
consider monojet + + photon topology in 13 TeV running with 300 fb“^. The curves show pro jections assuming a prompt 


photon. For BaBar projections, we propose a diphoton + ^ search (green dashed). See Sec. Ill for details. Also shown are 


constraints from a LEP diphoton + $ search |63| . LUX |64| . and Fermi line-searches |65| : see Sec. ^ We truncate LHC search 
results where the effective held theory dipole description is no longer valid, which accounts for sharp features at large dipole 
strength and large mass; see Sec. m for a detailed discussion. 


Inelastic Fermion Current 

Consider a familiar Dirac spinor t/j = iv charged 
under the U(l)£) gauge symmetry. The vector current is 

= T]'^a^r] — ( 6 ) 

where ry and ^ are two-component Weyl fermions. Gauge 
invariance allows only a Dirac mass term, rriD- However, 
when the symmetry is spontaneously broken, the compo¬ 
nents of Ip can also acquire Majorana masses: 

- + (7) 

Because the U(l)£) symmetry is restored in the limit that 
the Majorana masses go to zero, it is natural for the 


Majorana masses to be smaller than the Dirac mass. In 
this limit, the spectrum is split in the diagonal mass basis 
into two nearly equal Majorana mass eigenstates, where 
the mass eigenstates couple predominantly inelastically. 
For example, if the Majorana masses are equal, (m,( = 
mj), the mass eigenstates are xi = ~ O/v^ , X 2 = 

iv + with eigenvalues mi ^2 = tnn T wm, and the 

vector current now couples different states to one other, 

J'" = t{x\'^^X2 - xl^^Xi) = JiOM : ( 8 ) 

where the iDM subscript emphasizes the inelasticity of 
the interaction. 

In the more general case where niri ^ m^, the mass 
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eigenvalues are 


wi ,2 = \j~ 'm^y/A± {rrir, + m^)/2 , (9) 

and the physical splitting is A = mr^ + In this case, 
the vector current contains an additional elastic piece: 


= 


Tod 


+ 


m?jj + (toj — 

TOj — m 


m, 


)V4 


7'M 


vrj 


■y/4m|, + (mj - m^)2 


iDM 


(xJo-'"X2 - xlo-'^xi) ■ (10) 




Because we consider cases where A is not too much 
smaller than the the Dirac mass, we include the effects of 
the elastic coupling by assuming rrijj = A, = 0 when 
determining various constraints on the model. 


Xi 


-•,/\/\/\/\/\/\/\^ 7/.Z 


X2 


- 0 v/\/\/\/\/\/\/\^ ijz 



Inelastic Scalar Current 

FIG. 6: Leading self-annihilation (top) and inelastic annihi¬ 
lation (bottom) diagrams through the effective dipole inter- 
Similar conclusions follow if the DM is a complex scalar action depicted as a solid dot. 
ip with vector current 


= i[p*d^p-pd^p*). ( 11 ) 


B. y/Z-Mediated Dipole Interaction 


If the p scalar potential contains both [/(l)D-pi'eserving 
and -violating mass terms, the 1 mass terms for p can be 
written as 


- C D pL^P*p + ^P^PP + h-c.. 


( 12 ) 


where p and p are “Dirac” and “Majorana”-like mass 
terms for a scalar particle, which respectively pre¬ 
serve and break any gauge symmetry under which p is 
charged^. We have assumed that the p expectation value 
remains zero in order to preserve a symmetry to stabilize 
the DM. 


Diagonalizing the mass terms in Eq. (12 1 yields eigen¬ 


states = (</^ i p*)lxf^ with corresponding mass- 
squared eigenvalues ± The vector current is now 




(13) 


which is purely off-diagonal. Because the scalar mass 
eigenstates are real, no diagonal interactions with a single 
A' are allowed. The covariant derivative also yields diag¬ 
onal quartic terms of the form (jjitjiiA'A', which gives an 
elastic scattering mode involving pairs of gauge bosons, 
but the coupling to fermions in colliders and other exper¬ 
iments is suppressed by one loop. 


^ A mass term of type V(ip) = (p')^Im((p^)/2 is also allowed; how¬ 
ever, because the scalars are real in the mass basis, an SO(2) 
rotation among the scalars does not induce a diagonal coupling 
between A' and the mass eigenstates, and so it has no effect on 
our results. 


The other representative model we consider is mag¬ 
netic inelastic DM (MiDM) H5] . If we consider a DM 
Dirac fermion t/i = (ry ^1), it can have a direct coupling 
to the SM through the dipole interaction 

, (14) 

= + + , (15) 

where ps is a dipole moment of dimension inverse mass, 
a^^'' = and cr^(CT^) = (l,±crj), where ai is 

a Pauli matrix. This interaction is well-motivated and 
arises in a wide variety of DM models [iliniEDHTS]. 

In the absence of an exact DM U(l) symmetry, the 
components of ip can also have Majorana masses. In 
terms of the mass eigenstates, ^( 77 ) = (xi ±fx 2 )/'\/ 2 , this 
operator becomes 

+ h-c. , (16) 

where we have built 4-component Majorana fermions 
from each of the xi _2 Weyl spinors. For the dipole- 
mediated scenario, we assume this operator is the lowest- 
dimension interaction arising from a UV theory with 
additional (hyper) charged field content; for an explicit 
construction, see Ref. |40]. After electroweak symmetry 
breaking, this operator decomposes into 

^ = Y XWp.i"^"X 2 + ^XWp.^'^"X 2 + h.c. , (17) 

in the 7 and Z mass eigenbasis, where pzlp-y = 
— tan^vv; although higher-dimensional operators can 













Average Lab Frame Decay Length ^fs = 13 TeV 


modify this ratio |40j , we consider the simplest case where 
the ratio is given by the weak mixing angle, and we pa¬ 
rameterize our results in terms of Because dipole 
moments vanish for Majorana fermions, there is no di¬ 
agonal coupling akin to the interactions of fermion iDM 
with a dark photon. However, in computing DM anni¬ 
hilation rates, we include the effects of the higher-order 
t-channel processes XiXi 77; ZZ (see Fig.|^. 

As in the A! scenario, our focus is on covering the pa¬ 
rameter space that induces a thermal-relic annihilation 
in the early universe. However, a crucial difference in 
this case is that, for a given A, the annihilation rate 
depends only on and mi, so no additional assump¬ 
tions need to be made in order to compare different kinds 
of experimental bounds against the thermal relic bench¬ 
mark. Thermal Majorana MiDM with a mass splitting 
A > 0.2 mi is nearly excluded already by a combination 
of direct, indirect, and collider searches (see Fig. [^. For 
smaller mass splittings, only the mi > 100 GeV region 
is robustly ruled out by gamma-ray line searches, so the 
searches proposed in this paper are designed to target 
the remaining viable parameter space. 


III. COLLIDER SEARCH PROPOSALS 

Dark matter searches at high-energy colliders tradi¬ 
tionally feature missing (transverse) energy or MET) 
and fit into two broad categories: searches for DM pro¬ 
duced from the decays of additional new SM-charged 
states such as t-channel mediators [TTHIS], and searches 
targeting direct DM production through the reaction 
pp —>■ DM -I- DM -I- X where X is some visible SM state. 
The former class is more model-dependent by nature, al¬ 
though well-motivated frameworks like supersymmetry 
(SUSY) fall into this category; by contrast, the latter is 
more model-independent because it relies primarily on 
DM’s direct coupling to the SM. In recent years, there 
have been many proposed searches and analyses for DM 
pair-production, which yields in association with SM 
final states, including monojet, monophoton, and mono¬ 
boson P'81 - 139] . Indeed, the LHC is particularly well- 
suited for discovering classes of DM models with con¬ 
tact interactions of light DM, where the sensitivity of di¬ 
rect and indirect-detection experiments is sub-optimal, 
although these searches remain insensitive to contact- 
interaction strengths sufficient to induce thermal-relic an¬ 
nihilation rates in the early universe for light mediators. 

In an extended DS, two different particles in the DS 
can be produced in association at colliders, in contrast 
with mono-A searches that target only the production 
of ground-state particles. As a result, we show that the 
collider sensitivity to scenarios such as iDM can be en¬ 
hanced by tagging the decay products of the associated 
state(s), providing a powerful handle for background re¬ 
jection. We focus on the representative models in Sec.[Tl| 
However, there also exist models where the DM lives in 
an extended DS, which can give rise to more varied and 



FIG. 7: Lab frame decay lengths in the fermion iDM dark 
photon model for boosted X ‘2 de-exictation via xi 
which is mediated by an off-shell A!. Events are shown at 
= 13 TeV and requiring a leading jet px > 120 GeV. The 
results are normalized to j/ = 10“® for simple comparison with 
Figs.[2|-0 


spectacular signatures than those we consider (see, for 
example, Refs. [44H50] 1. 

We now propose a series of new searches at both the 
LHC and H-factories that can dramatically improve the 
sensitivity to the scenarios introduced in Sec. |H| We or¬ 
ganize our searches by model, since the models give very 
different signatures at colliders. 


A. Dark photon 

LHC 

We describe the LHC signatures of the dark photon 
iDM model introduced in Sec.|^ We focus on the regime 
of few-GeV iDM masses, where existing constraints on 
thermal DM are relatively weak and for which dedicated 
collider searches must be developed to tag the SM states 
from DM* decay. In our model, the excited state DM* 
decays via A'^*^ to DM + ff with the latter being SM 
fermions. Since Fdm* ~ (see Appendix]^, DM* 

is long-lived on collider scales for GeV-scale masses and 
moderate mass splittings (A/toi ~ O.OI — O.I), giving 
rise to decays within the LHC detectors at a displaced 
vertex. The signature is striking, but the leptons are 
typically both collimated and soft, motivating dedicated 
collider searches. We show a representative decay length 
distribution in Fig. A part of the parameter space 
was explored in Ref. [2l] in the context of iDM for fully 
hadronic X 2 Ai decays over a range of masses at hxed 
splitting, with a focus on contact operators with cou- 
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plings to quarks^. 

In our representative model, A' couples to the quarks 
and leptons via hypercharge kinetic mixing. We focus on 
the case where ttia' > m 2 + mi; in this scenario, the A! 
is produced on-shell in association with a jet^ and the 
A! decays to XiX 2 - The X 2 subsequently decays via an 
off-shell A! to x\ plus SM fermions. In particular, we 
focus on the decays to leptons, giving rise to a displaced 
dilepton vertex in association with a hard jet and missing 
momentum. In particular, we consider the reaction 

pp j + A'^ jx2Xi (18) 

ji^TxiXi- (19) 

The location of the displaced vertex determines the sen¬ 
sitivity of LHC searches to the A' model. If X 2 escapes 
the detector before decaying, the signature reverts to 
that probed by conventional monojet searches. If X 2 de¬ 
cays at distances < 1 mm, then large (quasi-)prompt SM 
electroweak, top, and QCD backgrounds become impor¬ 
tant. We find the best sensitivity is in between these two 
regimes. 

The signal also features distinctive kinematics for the 
dilepton pair. For mA’ ^ m^ ^ GeV, the A' is highly 
boosted when recoiling off the high-px jet. As a result, 
the decay products are highly collimated in the decays 
A' —)■ X 2 X 1 ^~^'(^~XiXi- III particular, the leptons are 
well within Ai? = -y/ Ar/^ + iXejA < 0.4 of one another 
(see Fig. I^for a representative benchmark point), giving 
rise to a lepton jet signature (see Refs. [TTl 1^[50ll80H86] 
for examples of lepton jets in other contexts). Similarly, 
the lepton jet will be collimated with the of tho event, 
as illustrated in Fig. In our analysis, we focus on the 
decays to muons since the backgrounds are smaller for 
this hnal state due to a lower fake rate. All of the above 
features of the signal combine to give a striking hnal state 
at the LHC, as we show below. 

The backgrounds for displaced vertex searches are 
small, but difficult to estimate due to contributions from 
rare heavy-flavor decays as well as random track crossings 
in the detector. Monte Carlo (MC) tools do not always 
model the backgrounds well, which is why the experimen¬ 
tal collaborations typically perform data-driven back¬ 
ground estimates for long-lived particle searches. While 
we lack the tools and the data to do such an estimate, in 
what follows we discuss in turn the most important SM 
processes that can mimic the signal and perform rudi¬ 
mentary estimates of the importance for each, arguing 
that our proposed hnal state could well result in a vir¬ 
tually background-free search. This conclusion is also 


® From this point on, we use X 2 and xi as placeholders for DM 
excited and ground states, respectively, regardless of whether 
DM is a scalar or fermion. 

^ We also include DM production via j + {Z ^ X2Xl)i although 
this is subdominant for the parameters we study. 
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FIG. 8: AR spectrum between the displaced dilepton tracks 
in the dark photon model. We consider A = O.lmi and A = 
0.4mi, assuming mi = 5 GeV and niA ~ 15 GeV. 

supported by the very low backgrounds observed by AT¬ 
LAS for soft displaced dilepton vertices [87]; the trigger 
and event selections for this analysis require at least one 
very high px lepton or trackless jet, however, and so the 
analysis in its current form is not sensitive to dark photon 
iDM signatures. The background estimates we provide 
must, of course, be verihed by detailed study by the ex¬ 
perimental collaborations. 

In our analysis, we perform all signal and background 
calculations in this section with Madgraph 5 |88|, with 
subsequent showering and hadronization in Pythia 6 
[55] . and jet clustering with Fast jet 3 [SO]- Moreover, 
we use Feynrules to define the two representative 
models in our study [5T] . 

Backgrounds: 

• Real photon conversion. Non-prompt leptons can 
arise from the conversion of real photons via colli¬ 
sions with material or gas in the detector. A photon 
traversing the inner detector (ID) inside ATLAS or 
CMS will go through one radiation length of ma¬ 
terial, so an O(l)-fraction of the photons convert 
into an e“'"e“ pair. Moreover, the production cross 
section for pp —>• jx{Z vv) at = 13 TeV is 

« 100 fb, after requiring a jet with px > 120 GeV 
and a sufficiently energetic, isolated photon. Nev¬ 
ertheless, there are several handles to significantly 
reduce this background. First, events where the 
photon carries some of the momentum away can 
be rejected by applying stringent isolation require¬ 
ments on the leptons. Second, leptons originating 
from material-dense parts of the detector and/or 
consistent with photon conversion can be vetoed. 
Third, the invariant mass distribution of the lep- 
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FIG. 9: Distributions of | A0| between the momentum of the lepton jet and the of the event in the dark photon model. We 
consider A = O.lmi and A = 0.4mi, assuming mi = 5 GeV and m^/ = 15 GeV. 


tons from photon conversions will peak at zero mass 

as they originate from an on-shell photon. Finally, 

the photon conversion probability to muons is sup- 

2 

pressed relative to electrons by ^ « 10“^. All of 

these considerations combine to allow an estimate 
of negligible background from photon conversions. 

• QCD. Displaced tracks could originate from QCD- 
initiated jets, particularly those giving rise to long- 
lived B or K hadrons which in turn decay to tt 
and/or /i. Estimating this background from first 
principles is not feasible due to the dependence 
on hadronization effects and the challenge of esti¬ 
mating muon mis-identification rates; nonetheless, 
we determine an approximate upper bound on the 
probability for a QCD-initiated event to give a hard 
leading jet with px > 120 GeV, and two muon- 
tracks with pt > 5 GeV appearing at a displaced 
vertex (both muon tracks have transverse impact 
parameter |fio| between 1 mm and 30 cm, and the 
point of closest approach of the tracks is < 1 mm). 
We estimate this probability to be < 10“^, which 
bounds the QGD cross section to be < 100 fb. This 
requirement is before requiring significant missing 
energy from hadrons in the event, and before re¬ 
quiring that the missing energy be near the muon- 
tracks or any other kinematic features characteris¬ 
tic of signal. 

• Pile-up. In the QCD estimate, we assumed that 
the jet and the muons from long-lived hadron de¬ 
cays originate from the same primary vertex. In 
upcoming running of the LHC, the increased lu¬ 
minosity comes with the price of a large number 
of primary vertices per bunch crossing due to pile- 


up. Therefore, it is possible that the soft, displaced 
muons could originate from a different primary ver¬ 
tex. Since the signal muons are highly collimated, 
however, they point in the same direction as the 
long-lived particle, which passes through the pri¬ 
mary vertex with the high-px jet- Therefore, even 
though each muon has a high impact parameter, 
applying a selection requirement that the dimuon 
momentum approximately point back to the pri¬ 
mary vertex can be used to suppress long-lived 
hadronic backgrounds from other primary vertices. 

• Jet + di-tau. The cross section for a high-px jet, 
along with two r leptons within Ai? < 0.01 of one 
another, at the LHC at -y/s = 13 TeV is « 10 fb. 
Accounting for the requirement that both taus de¬ 
cay to a muon further reduces this rate to ~ 10“^ 
fb. In addition, for this background component to 
mimic the signal, both taus need to decay within 
~ 100 pm of each other. And finally, we note that 
since each muon-track originates from a different r 
parent, the distribution will be distinct from 
the signal where both tracks originate from the X 2 ■ 

• Jets V —>■ A. potential background may 

originate from the reaction pp —jets -I- V, with 
V either a Z or a tV boson decaying to give miss¬ 
ing energy. For this background to contaminate 
the signal region, one would need the two tracks to 
originate from the jets. Through a reasoning anal¬ 
ogous to that used above for the QCD background, 
this background component should be in the range 
of less than ^ 0.01 fb, and so relatively negligible 
for our analyses. 

• Backgrounds from fake missing energy. Typically, 
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experimental analyses require a minimum separa¬ 
tion between missing energy and other objects in 
the event to suppress fake missing energy from 
calorimeter or momentum mis-measurement. By 
contrast, our signal is collimated with the miss¬ 
ing energy, and so fake missing energy is a poten¬ 
tial concern. We exploit the fact that the muons 
from signal decays are relatively soft (typically, the 
summed muon pt/.^t is 0.2; see Fig. 10), and 
so fake missing momentum is not expected to be 
important. 


Signal region: The above considerations motivate the 
following selections for the signal region: 

• Trigger on a mono jet -I- For Run 1, for in¬ 

stance, CMS used a > 120 GeV trigger [5^ . 
where th® missing momentum as computed 

in all subsystems excluding the muon system. We 
assume such a trigger for our sensitivity estimates 
for y/s = 13 TeV, although note that the exact val¬ 
ues for Run 2 could be higher. The additional use of 
the soft leptons could help keep trigger thresholds 
low; for example, ATLAS has an analysis which 
has requirements as low as px > 6 GeV for muons 
at trigger level (92]. Nevertheless, we also checked 
that with a trigger of ptj > 200 GeV andA^r > 200 
GeV, the signal sensitivity is degraded by approxi¬ 
mately a factor of two in rate; 

• One leading jet with pT > 120 GeV and allow only 
one extra jet with px > 30 GeV; the leading jet 
and ^x should be back-to-back; 

• One displaced muon jet, pj, consisting of at least 
two muons with |do| between 1 mm and 30 cm, 
and whose tracks cross within 1 mm; the two muon 
tracks each have px > 5 GeV; 

• $T> 120 GeV; 

• <0.4. 

We show our projections for the LHC sensitivity to this 
topology at -/s = 13 TeV, assuming C = 300 fb“^ of in¬ 
tegrated luminosity in Fig. [^through Fig.|^for fermionic 
and scalar DM. As motivated above, we assume a mostly 
background-free signal region, and therefore plot sensi¬ 
tivities for 10 signal events. The unique kinematics of 
the signals we study could allow the LHC to probe the 
very couplings responsible for establishing the DM abun¬ 
dance in the early universe through thermal freeze-out. 


B factories 

H-factories have the potential to make significant 
progress in the exploration of DM with inelastic inter¬ 
actions for DM masses within kinematic reach (below a 


few GeV). This is due to the very high luminosity and 
clean environment of an collider. 

There are two potential avenues to pursue at a B- 
factory. One possibility is to search for direct production 
of the ground and excited states with subsequent decay 
of the excited state to (displaced) tracks. It is unclear, 
however, whether such displaced tracks are sufficiently 
energetic or well-reconstructed to pass the trigger in 
analyses such as Ref. [33]. Alternatively, one can trigger 
on and reconstruct a visible SM state Y produced 
in association with DM (e<e“ —>■ DM -|- DM* -|- Y). 
In particular, the BaBar experiment instrumented a 
monophoton trigger for « 60 fb“^ of the total dataset. 
In our study, we conservatively consider only the latter 
scenario as a trigger for our proposals, although both 
possibilities should be investigated by H-factories. We 
base the following results on the analysis from Ref. [94j . 
which used a photon trigger with threshold > 2 GeV. 

Monophoton + missing energy: The analysis from 
Ref. |33| performed a search for (untagged) decays of 
T(3S') —A° -I- 7 (£ « 25 fb“^), where A° is an invisibly- 
decaying pseudoscalar, with a stringent veto on addi¬ 
tional activity in the detector. The dark photon in our 
model is produced through the reaction e+e“ —>■ 7 -I- A', 
with subsequent decay A' —)■ X2Xi- Although our model 
produces visible states in X 2 decays, the kinematics of the 
dark photon signal can still populate the BaBar signal 
region, which consisted of a bump search in the missing 
mass variable ~ 2A.y_cMWx(3S)- For our 

signal to appear in this search, the X 2 has to decay either 
outside the detector or into soft final states that fall be¬ 
low BaBar’s thresholds (we use the thresholds listed in 
Ref. [95]). 

A complication of this analysis is that, for tua' < 
1 GeV, the signal could appear in the BaBar control 
region, in which case the signal mimics the kinematics 
of the irreducible 77 background. In this mass regime, 
we set a conservative bound by assuming that the signal 
represents all of the events in the control region. 

The results of the BaBar monophoton -|- missing 
energy recast are shown in Fig. [^ to [^ We also provide 
a projection for Belle H assuming 50 ab“^ of luminosity 
with an instrumented monophoton trigger [96j . For a 
more extensive discussion of the details of this analysis, 
see Ref. m and Ref. [33] . 

Monophoton + displaced tracks + missing en¬ 
ergy: A potentially more striking signature of iDM at 
H-factories could be uncovered by a re-analysis of BaBar 
data. In particular, the reaction > 7 -I-A' can give 

rise to displaced tracks and missing momentum in the 
final state. As before, we assume use of the monophoton 
trigger, and offline selection of two displaced leptons with 
p > 100 MeV and transverse impact parameter |do| be¬ 
tween 1 cm and 50 cm, as in Ref. [33] . Based on Ref. [98] , 
we use a lepton reconstruction efficiency of 50% in our 
estimates. 
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FIG. 10: pt distributions for the leading and subleading muons in X 2 —>■ decays at the LHC at y/s = 13 TeV in the 

A'-mediated scenario for representative masses and splittings. 


We expect backgrounds from resonances (arising from 
the decay of a hadron, or through radiative return) 
in this channel to be low, particularly after requiring 
significant and removing dilepton pairs consistent with 
hadronic resonances. Studies of backgrounds in related 
searches for —)• J /^/>7 [S5] and radiative decays of 

T —> "f{AP —)• /t+zr”) | 100 j (where A° is a light exotic 
scalar) suggest that this may indeed be the case for our 
proposed channel. Additionally, from Ref. m. another 
potential background is that from 77 r"'' 7 r“, and 77 with 
one of the 7 converting to a pair. However, the 

former can be reduced with the requirement that the 
tracks originate from a high impact parameter vertex, 
and the both the former and the latter could be reduced 
through a combination of a missing mass cut and a cut 
on the invariant mass of the tracks. 


Results: The above proposed searches at BaBar prove 
complementary to the searches at the LHC that we ad¬ 
vocate. In particular, we find they have the poten¬ 
tial to cover thermal-relic territory for the 0 ( 10 )% frac¬ 
tional mass splittings that are the focus of our analysis. 
Figs. i-i illustrate the potentially powerful reach that 
BaBar could achieve with a dedicated monophoton -|- dis¬ 
placed tracks search. Additional improvements could be 
achieved by future B-factories |96j depending on whether 
or not they are instrumented with a monophoton trigger, 
especially outside of the control region where the sensitiv¬ 
ity scales as >/£; therefore, our analysis provides further 
motivation for the development of a monophoton trigger 
for Belle B. 


B. Magnetic Dipole Interaction 

LHC 

The second simplified model we consider is dark matter 
coupled inelastically via a magnetic dipole moment (see 
Sec. E- In this scenario, the excited DM state X 2 decays 
via X 2 —>■ Xi + 7 - We are interested specifically in the 
^ 100 MeV-100 GeV, 0(10%) splitting inelastic limit 
considered earlier. As before, the decay products of X 2 —> 
Xi -b 7 are typically too soft to serve as the main trigger 
objects, and so we rely on the associated production of 
a high-px jet. Thus, we predict a pp ^ j + + 7 

signature. Existing work has studied the scenario with 
a hard photon originating from larger splittings between 
DM states in both the prompt and long-lived limits [iQl 
EZ]. 

There are two principal distinctions between the dipole 
scenario and the dark photon considered earlier. The 
dipole is a dimension-5 operator, and so the decay width 
of X 2 through the dipole in the limit of small splittings 
A goes like T ^ (see Appendix [A|); by contrast, de¬ 

cays through an off-shell dark photon scale like jm\, 
and is suppressed by 3-body phase space. As a result, 
the decays are prompt over a wide range of the dipole 
parameter space, and consequently the backgrounds are 
significantly larger than in the displaced muon jet anal¬ 
ysis. Furthermore, it is more challenging to reconstruct 
soft photons than soft muons, with the photon recon¬ 
struction efficiency > 0.5 only above Et = 15 GeV (see, 
for example. Ref. una). Thus, the sensitivity of a dedi¬ 
cated search for the existence and kinematics of the soft 
photon is lower than for the dimuons. Nevertheless, we 
find that dedicated monojet -I- photon -|- missing energy 
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searches can be competitive with existing limits and cover 
interesting, unexplored parts of the parameter space, par¬ 
ticularly for i3-factory probes. 

We consider two cases: in the first, X 2 decays on 
lengths shorter than the pointing resolution of the elec¬ 
tromagnetic calorimeter (ECAL), namely < 25 mm for 
E'Y-y ~ 20 GeV |103j . This region of lifetimes can be 
probed by dressing the canonical monojet searches with 
a prompt, soft photon requirement. Such searches have 
significant electroweak backgrounds, but our proposed 
searches are competitive with existing constraints and 
can have some sensitivity over new parameter space due 
to the distinctive signal kinematics. Second, for small 
dipole moments, the decay occurs on macroscopic scales 
and the photon from X 2 decay does not point back to 
the primary vertex. Such non-pointing photon signatures 
have been considered in SUSY searches m for pairs of 
photons. For our signal, X 2 is highly boosted if the pho¬ 
ton is hard enough to be seen in the detector, and so it 
still points approximately in the direction of the primary 
vertex; consequently, when the photon decays before the 
ECAL, its pointing displacement is typically ^ cm, and 
so cannot be effectively distinguished from prompt de¬ 
cays. Thus, we focus on prompt searches here. 

We begin by examining the relevant backgrounds. 
In our analysis, we use the same MC tools as for the 
dark photon model described above. We re-normalize 
backgrounds using an average next-to-leading order 
AT-factor calculated using MCFM |1041 1105) : as we show 
below, the signal over most parts of parameter space 
is produced via {Z —)■ X 2 X 1 ) + j, and so we apply the 
Z-l-jets AT-factor to signal as well. 

Backgrounds: 

• Jets + "f + £ + 12 . This dominant background 

comes from kT- 1 -jets, where the lepton undergoes 
bremsstrahlung and emits a photon. While lepton 
vetoes can be used to suppress this background, a 
hard photon and soft lepton occur in a substan¬ 
tial fraction of events. The resulting photon is cor¬ 
related with the direction of the lepton, and for 
boosted W bosons, this results in a soft photon 
aligned with This provides an irreducible back¬ 
ground mimicking our signal. We assume that pho¬ 
tons produced by collisions of the lepton with gas 
and detector material can be suppressed through 
careful event selection, and so the dominant back¬ 
ground comes from prompt photons. Applying the 
lepton veto from the ATLAS analysis [3H] and a 
sample of preselection requirements {pTj > 120 
GeV, > 100 GeV, an isolated > 15 GeV), 
and 7 )] < 1.5, this background has a cross 

section of approximately 250 fb at ^/s = 13 TeV. 

• Jet + (Z ^ Dv). This is the dominant back¬ 
ground for monojet searches at the LHC, and it pro¬ 
duces photons through decays of pions and other 
hadrons. However, the photons are typically nei¬ 


ther energetic nor isolated, and the presence of 
nearby hadronic activity can be used to greatly sup¬ 
press this background. After requiring the same 
preselection cuts as for W-l-jets, the cross section is 
approximately 5 fb. The photon is typically aligned 
with the jet, and unlike for the signal and W + jets 
background, this can be used to provide significant 
discrimination. There may also be fakes where an 
additional jet fakes a photon, although we expect 
that they can be suppressed with sufficiently tight 
photon identification requirements so that the rate 
is subdominant to the W-l-jets background. 

• Jet + x+iZ —)■ Piy). This background is the same as 
the above but with the inclusion of a hard, prompt 
photon. Its cross section after pre-selection is ap¬ 
proximately 50 fb, which is larger than the pure 
.Z-l-jets rate. The Z-|-jets background has very sim¬ 
ilar kinematics to the .Z-l- 7 -l-jets, with the photon 
typically correlated with the direction of the jet. 

• Fake missing energy. Because the photon is aligned 
with the missing energy, it is possible that QCD 
backgrounds fake both missing energy from jet mis- 
measurement as well as a photon from the same jet. 
We cannot simulate this background. However, the 
ratio of E-Y--jlt‘T is typically small (< 20 %) for sig¬ 
nal and large for fake backgrounds, and so we ex¬ 
pect that it can be greatly suppressed through 
isolation requirements. We expect that this back¬ 
ground should be negligible, and can be suppressed 
by increasing the px and thresholds if necessary. 

While there are some fake backgrounds that we 
cannot simulate, the above estimates are sufficient for a 
qualitative estimate of the parameter space reach of a 
monojet -|- photon analysis. 

Signal region: We employ a monojet -I- missing en¬ 
ergy trigger. At 8 TeV, the lowest threshold mono¬ 
jet -I- .^x trigger required monojet px > 80 GeV and 
4^x > 105 GeV [3S] , with the efficiency plateau somewhat 
higher. This threshold may increase at higher energies; 
however, if the photon from X 2 decay is sufficiently hard, 
it may be used to lower the mono jet trigger threshold®. 
Therefore, we perform an analysis with this trigger as the 
minimum threshold. We also briefly consider the effect 
of a higher threshold: for example, if we consider instead 
a minimum threshold of pxj > 200 GeV and .^x > 200 
GeV, the results are nearly unchanged for A < O.lmi 
because, in this regime, the X 2 needs a large boost for 
the photon to mass minimum reconstruction thresholds, 
driving the signal events to the large px and ^-x regime. 


® In the CMS high-level trigger at 8 TeV, photons as soft as Ex > 
23 GeV were used to keep the diphoton trigger threshold low 

[irob . 
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For larger A, less boost is needed for signal events to pass 
the photon cuts, and we find a degradation of < 25% for 
the prompt search. This does not qualitatively change 
our results. 

The backgrounds presented above limit the effective¬ 
ness of the prompt analysis and motivate the signal re¬ 
gion for our analysis, which is obtained by cutting on 
various kinematic observables. We scan over values of 
the cuts, and the least constraining values of each cut 
are listed below: 

• One central hard jet (|r?| < 2.5, px > 120 GeV) 
and allow only one extra jet with px > 30 GeV; 
the leading jet and^x should be back-to-back; 

• Missing momentum > 105 GeV; 

• One isolated® central photon (minimum ifx > 
15 GeV). We apply a flat 50% identification ef¬ 
ficiency for the photon, which is appropriate for 
Sx « 15 GeV 

• Difference in azimuthal angle |A(/)(.^X)7)1 < 1) as 
motivated by Fig. El 

• Transve rse mass: _ 

mx = — cos A^(.^X) 7 )] < 50 GeV; 

• The long-lived particle decay must occur before 
reaching the EGAL and faster than 4 ns. We also 
impose a pointing requirement: we calculate the 
distance |Az.^| by using the photon direction to 
determine the point of closest approach to the pri¬ 
mary vertex in the transverse plane, and then defin¬ 
ing |Az-y| to be the longitudinal distance between 
this point and the primary vertex |103j . We require 
that |Az.y| be within the detector resolution, and so 
satisfies \^z^\ < 15 mm 1^/50 GeV/ifx^ . 

The transverse mass and At); observables are some¬ 
what correlated, but we find performance is up to 25% 
better with both included. For each model benchmark 
point, we optimize over simple cuts in each observable, 
and use this to determine the 2cr reach at 300 fb“^ 
assuming a systematic uncertainty of 10% (and requiring 
a minimum number of three signal events). 

Results: We show our results in Fig. under “LHG 
prompt”. Over the range of splittings considered, our 
proposed search is competitive with existing bounds 
and can constrain open parameter space for MiDM. The 
cutoff of the sensitivity at mi ^ 50 GeV is not due to 
any kinematic effect, but rather due to a breakdown 
of the MiDM effective theory as we discuss below. 
Nevertheless, we expect the LHG to have reach for 


We require no more than 4 GeV of energy deposited within AR < 
0.4 of the photon, not including the photon itself. 


higher masses in a UV-complete theory. 

Validity of the Dipole Effective Field Theory: The 

dipole operator is simply the lowest term in a series ex¬ 
pansion generated by a UV complete theory with new 
mass states at some scale A. When the momentum flow¬ 
ing through the photon, , becomes larger than the scale 
of the new mass states, then the new states can appear 
on-shell and invalidate the effective field theory (EFT). 
In the case of the dipole operator, the largest scale A 
at which new states appear is A « g/ fi-y, which corre¬ 
sponds to UV physics that saturates the perturbativity 
limit. Therefore, our analysis necessarily breaks down 
near ^ A^, and processes with > A^ are not valid 
unless resolved by a form factor. 

To avoid sensitivity to an unphysical regime, we veto 
events where the dipole EFT description is invalid in 
calculating our results shown in Fig. [38l 11071 I108j . 
The valid MG events in our analysis often come from 
[Z —>■ X 2 X 1 ) +U which explains the cutoff in sensitiv¬ 
ity due to EFT invalidity for mi > mzj^. in Fig. A 
study of Dirac MiDM models in a particular UV comple¬ 
tion m found that the EFT description was conservative 
for rates with px lower than the cutoff, since it neglects 
the typically large contribution to the cross section for 

^ (this also persists at two loops |109j V. Since our 
proposed analysis has excellent sensitivity even in this 
conservative limit, and the regions of EFT invalidity are 
often well-covered by other constraints from Sec. |V] we 
present our results in the EFT limit. It is still an inter¬ 
esting problem to understand the reach of our search in 
a particular UV complete theory. 

Additionally, in the limit where the EFT is not always 
valid at the scale of LHC parton-level interactions, direct 
searches for electrically charged states that induce Rj can 
become important [nnum], but charged states below 
100 GeV are severely constrained by LEP. There may 
be overlap between such direct searches and the MiDM 
constraints for large dipoles, but the constraints from 
charged state searches are model-dependent. 


B factories 

As with the dark photon, the H-factories can have 
impressive coverage of the MiDM thermal relic param¬ 
eter space for DM masses below a few GeV. If the 
photon from the X 2 decay registers in the detector, then 


^ For ^ A^, the EFT description is far more optimistic than 
the UV complete theory, which in part explains our decision to 
employ a method of truncation. Furthermore, the authors of 
Ref. 1471 found that, even with the large discrepancies in the 
spectra between the EFT and UV completion, the difference in 
coupling sensitivity of the different methods is 0(25 — 50%). 
Given the large uncertainties in other aspects of our proposed 
analysis, this is within the tolerance of our projections. 
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MiDM Backgrounds MiDM Signals, A=0.1 mi 




FIG. 11: Distributions of |A(;A| between the momentum of the photon and the of the event for SM backgrounds (left) and 
signal in the dipole model (right). For signal, we consider A = O.lmi, mi = 5 GeV and 50 GeV. 


signatures of MiDM can appear in prompt or displaced 
photon searches via e+e” —7 + X 2 + Xi- before, we 
conservatively assume that events are recorded via the 
BaBar monophoton trigger. 

Two prompt photons and missing energy: For 

decay lengths < 1 cm, the photon from the decay 
X 2 Xi + 7 appears as a prompt photon. Thus in this 
case, the signal of interest is two prompt photons, one 
of which is used to pass the monophoton trigger, and 
missing energy arising from the decay of the neutral X 2 - 
In particular, we consider a signal region with a leading 
photon satisfying > 2 GeV to pass the trigger, a sub¬ 
leading photon with E > 20 MeV and missing energy 
.^ > 50 MeV. This scenario inherits low backgrounds 
from the monophoton search, and we assume that such 
a signal region is background-free in our projections in 
Fig.H as the dominant background is 77 and can be 
reduced with a missing momentum cut. Thus, we use 10 
signal events as a benchmark for the sensitivity to this 
model. Note that we apply the same acceptance cuts on 
the leading photon —0.31 < cos 6* < 0.6 (calculated in 
the CM frame) as Ref. [M]. This centrality requirement 
was important for the analysis from Ref. |94j . since the 
energy resolution is degraded at smaller 9*, but less 
crucial for the MiDM signal as this final state no longer 
calls for a bump search. Relaxing this requirement could 
indeed result in a stronger sensitivity to this scenario. 

One prompt photon, one non-pointing photon 
and missing energy: An orthogonal final state to con¬ 
sider has a decay length > 1 cm. We conservatively 
only consider scenarios where the decay occurs within 
£^2 < 50 cm. The signal region we define in this sce¬ 
nario is one prompt photon with E^ > 2 GeV, and one 


non-pointing photon with E^ > 20 MeV. As for the case 
above, we assume a background-free scenario in this sig¬ 
nal region for our projection; we find the signal sensitivity 
to be sub-dominant to the diphoton -|- $ signal, although 
it may be useful in the case of unexpected backgrounds 
for the prompt search. 


IV. INELASTIC DM RELIC ABUNDANCES 

In the standard thermal freeze-out scenario, the DM 
abundance today is fixed by the DM number density at 
the time that its annihilation rate into SM states falls 
below the Hubble expansion rate. For any given model, 
the observed relic abundance singles out a particular re¬ 
gion in parameter space that is an obvious target for ex¬ 
perimental searches. In many models, the thermal relic 
parameters are challenging to directly probe at colliders. 

As we have seen, inelastic DM models offer additional 
handles at colliders to probe small couplings. To assess 
the sensitivity to the thermal relic cross section, we dis¬ 
cuss in this section the calculation of the DM abundance 
in the (M)iDM models from Sec. |nj leading to the con¬ 
clusion that thermal scenarios can be probed in many of 
the searches proposed in Sec. |III| In iDM models, DM 
annihilation proceeds predominantly via co-annihilation 
of the mass eigenstates; since X 2 is heavier, its abundance 
depletes faster than xi due to decays (x 2 —>■ Xi + SM) 
and scattering (x 2 + SM —Xi + SM), and so freeze-out 
occurs earlier for mass splittings larger than the freeze- 
out temperature. This suggests that thermal relic iDM 
models feature larger couplings to compensate for the less 
efficient freeze-out. 

In terms of the variable x = TO 2 /T, the coupled Boltz¬ 
mann equations for the DM densities V (normalized to 
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the entropy density) are: 
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( 20 ) 


where Yi = ni /s is the comoving number density of each 
species, a ( 0 ) superscript denotes an equilibrium quan¬ 
tity, s(r) = 27r^gs_*T^/45 is the entropy density, and 
Aa, As, and 7/5 are dimensionless annihilation, scatter¬ 
ing, and decay rates respectively. gs,*(T) is the number 
of entropic degrees of freedom. The first line of the right- 
hand side characterizes the change in DM density due to 
co-annihilation, the second line gives the change due to 
self-annihilation, and the third and fourth lines charac¬ 
terize scattering and decay processes that keep xi and X 2 
in chemical equilibrium with one another and in kinetic 
equilibrium with the SM. Using the Hubble rate during 
radiation domination H{T) = 1.66^^T^/mp£ (g* is the 
number of relativistic degrees of freedom), the dimen¬ 
sionless rates are defined to be 




( 21 ) 

As 

- Xlf)) 

( 22 ) 

Id 

= H{m2)~^{T{x2-)■ X1+SM)), 

(23) 


for X 1 X 2 co-annihilation, X 2 f Xif inelastic scattering, 
and X 2 —>■ Xi + SM decays respectively. The diagonal 
rate A^ is non-zero if there exist processes that allow 
XiXi SM -|- SM annihilation. 

For the dark photon model, the scalar dark matter 
scenario is purely inelastic and so A^ = 0. For fermion 
DM, there exists a self-annihilation channel whose rate 
is proportional to the difference of Majorana masses in 


Eq. (10), and is also p-wave (helicity) suppressed for the 


SM vector (SM axial) current. For the pure dipole sce¬ 
nario, the XiXi 77; jY, and ZZ channels are always 
open if kinematically accessible, but the self-annihilation 
rate is suppressed by additional powers of the dipole mo¬ 
ment. 

As in most co-annihilation scenarios, the scatter¬ 
ing/decay processes preserve kinetic and chemical equi¬ 
librium between X 2 and xi throughout freeze-out, and 
so the system of Boltzmann equations for Yi 2 can be re¬ 
placed by a single Boltzmann equation for Ytot = Yi+Y 2 , 


dYto, 

dx 


y 2 
^ tot 


L ( yW + F 2''^^)2 


70)^ 


- 1 


E$Yr>rr- («) 


* j 


iDM Thermal Freeze-Out 



X = m2jT 


FIG. 12: Freeze-out for fermion iDM (including co¬ 
annihilation and sub-dominant self-annihilation) mediated by 
an s-channel A' with m\ = 10 GeV, A = 0.2mi, and 
mA' = 3 mi with {av) ~ 10~^^cm^ s“^, for which ~ Ddm 
at late times. The solid (dashed) curves represent the actnal 
(equilibrium) number densities for the xi ,2 species and we de¬ 
fine the dimensionless evolution parameter x = m^jT. Note 
that the excited state continues to steadily decay and down- 
scatter into Xi off SM particles even after xi has frozen out. 


This approximation is valid over our parameter space. 

Considering an example point in the dark photon 
model, we show in Fig. the xi and X 2 yields as a 
function of m 2 /T. For each model, we determine the pa¬ 
rameters of the theory that give the observed DM relic 
abundance as a function of mi, and we show these curves 
in Figs.[2][^ We provide more comprehensive information 
on the rates that appear in the Boltzmann equations in 
Appendix 


V. CONSTRAINTS 

In this section, we consider other constraints on the 
parameter space of the dark photon and dipole mod¬ 
els, reviewing those which are complementary to collider 
searches and those which are ineffective in iDM models. 
These probes include direct detection experiments, preci¬ 
sion measurements of SM parameters, indirect detection, 
and LEP. 


A. Precision Electroweak and QED Measurements 

For models with new neutral gauge interactions, 
mixing between the massive gauge bosons can lead to 
shifts in observed SM electroweak couplings that are 
excluded by electroweak precision and other observables. 
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Dark photon: At high energies, the gauge invariant ki¬ 
netic mixing term is between the A and hypercharge. At 
low energies, this leads to mixing with the Z and can shift 
the observed mass at the Z pole to an unacceptable ex¬ 
tent. Precision electroweak measurements are therefore 
sensitive to A independent of its decay modes. Con¬ 
straints from LEP [62j impose a generic constraint on 
e < 10“^ at low mA' and a stronger limit near the Z 
mass. We show the exclusion region in Figs. 

There are additional constraints from the observed 
value of the muon’s anomalous magnetic dipole moment, 
{g — 2 )^. Despite the ^ 3cr discrepancy between theory 
and observation of {g — 2 )^, the ba bounds can still be 
used to set a conservative limit on the e/m^' ratio. In 
terms of = (5 — 2 )^/ 2 , a kinetically mixed gauge boson 
induces a shift [5] 


Attfj, 


ae^ ^ 2z(l — z)^ 

27 ^ Jo (1 - z)^ + {niA'/mfif z 


(25) 


which constrains the upper-left portion of the parameter 
space in Figs [^“1^ 


The LEP searches have also been used to constrain mod¬ 
els with dipole transitions between multiple Dirac DM 
states, giving rise to hard photons m- 


C. Indirect Detection 


Null results in searches for DM annihilation in cosmic 
rays or at the epoch of recombination are typically very 
constraining of elastic thermal DM models with s-wave 
annihilation to states other than neutrinos; for example, 
measurements of the cosmic microwave background 
(CMB) typically exclude such models for DM masses 
below ~ 10 GeV |11214114] . In iDM models, however, the 
relic abundance is largely set by co-annihilation of DM 
states in the early universe, whereas the X 2 abundance 
is depleted today, suppressing indirect detection signals; 
the mass splittings we consider are too large for signals 
from collisional excitation m- Nevertheless, sub¬ 
leading processes may mediate residual xi annihilation 
today, and we examine the constraints. 


MiDM: Since MiDM does not introduce any new gauge 
boson eigenstates, there are no relevant constraints on 
the coupling from precision SM measurements. 


B. LEP Constraints 

Dark photon: iDM can be produced via Z — >■ X 2 X 1 , 
as well as via radiative return production of A —>■ X 2 Xi- 
We find that the LEP searches for soft leptons and miss¬ 
ing energy are subdominant to the Z-mass constraint 
discussed earlier. 

MiDM: We considered several constraints, including the 
Z invisible width and direct searches for final states with 
photons. Over much of the parameter space, we find 
that the most powerful constraint comes from a LEP-1 
search for two photons -I- missing energy [63]. There is 
excellent sensitivity to MiDM due to the relatively low 
photon thresholds (two central photons with > 1) 
and the very low backgrounds (no events observed with 
rriyj > 5 GeV). We also include the constraint arising 
from the Z —>■ X 2 X 1 contribution to the total width [Tj, 
which is particularly important for fractional mass split¬ 
tings below 10% where the diphoton -|- MET search is 
less effective® . 

The excluded region is shown in Fig. [^ covering parts 
of the thermal relic line for larger splittings (A > 0.2 mi). 


The photon from X 2 ^ Xi + 7 decay is often not soft enough 
to avoid the veto used in Z invisible decay searches, and so Z 
invisible width constraints are typically subdominant to the total 
width measurement. 


Dark photon: In the scenarios we consider, mi < m^', 
and in the limit of zero splitting, self-annihilation only 
occurs in either four-body final states or loop-induced 
processes, both of which lead to tiny indirect detection 
cross sections. If A is larger, self-annihilation can also 
occur through p-wave-suppressed (via mixing with 7 or 
Z) or helicity-suppressed (via mixing with Z) operators, 
both of which give suppressions that render indirect 
detection signals negligible. Therefore, there are no 
relevant constraints on this scenario. 

MiDM: Processes with two insertions of the dipole oper¬ 
ator allow for the tree-level annihilation XiXi 11 
xZ (see Fig. 1^ [40]. Dark matter annihilating via this 
process gives rise to mono-energetic gamma rays origi¬ 
nating from the Galactic Center. Even though the indi¬ 
rect detection cross section is suppressed relative to the 
single-photon process determining the relic abundance, 
gamma ray line searches have sensitivity to sub-thermal 
cross sections over a broad range of masses. We recast 
the limits from Fermi’s line search. Ref. |65j . and show 
the result in Fig.[^ There is considerable uncertainty in 
the DM annihilation rate due to the unknown halo pro¬ 
file. To avoid over-stating the strength of the indirect 
detection bounds relative to collider searches, we use the 
Isothermal profile bounds in our recast, which is a factor 
of few weaker than an NEW profile. 

Measurements of the CMB also place constraints on 
DM annihilation. DM annihilations near redshift z ~ 10® 
can inject energy into the visible sector and ionize hy¬ 
drogen at the CMB’s surface of last scattering HIMIi]. 
However, we find that current constraints from photon 
line searches are stronger than those from the CMB, and 
so we do not include CMB constraints in our plots. 










D. LHC 
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Searches during Run 1 at the LHC in principle can 
already set constraints on the representative models of 
iDM we consider. Here we discuss them in turn. 

Dark photon: We checked that standard 

monojet+X+i7T searches for DM at the LHC don’t 
constrain new parameter space for the dark photon 
iDM model. However, a recent CMS analysis aimed at 
stop pair production in SUSY which searched for low 
multiplicity of jets, in addition to a pair of soft prompt 
dileptons and missing energy, sets limits on the param¬ 
eter space of our model |116j . In particular we find it 
provides complementary coverage to the searches we pro¬ 
pose, namely for larger mi and A, as shown in Figs. 

MiDM: We have verified that potentially constraining 
searches such as a diphoton and missing energy analysis 
by ATLAS from Ref. |117j do not cover new parameter 
space in the MiDM model. 


Xi Xi 



FIG. 13: Leading elastic scattering diagram for DM that 
couples off-diagonally to a vector mediator or to a photon 
via the dipole interaction in the MiDM scenario. For A much 
larger than the kinetic energy, we can contract the X 2 internal 
line and use the effective interaction shown in Fig. |15| 


E. Direct Detection 


One of the principal motivations for iDM models is 
that xi-nucleus tree-level elastic scattering XiA —)• XiA 
is no longer the dominant scattering reaction; the corre¬ 
sponding inelastic process, Xi-^ X 2 N is only allowed 
for certain kinematic configurations, modifying the naiVe 
predictions of DM spin-independent scattering [T2]. In 
particular, when A 3> |(f p/2Mjv for momentum transfer 
Iql and target mass Mjv, the inelastic process is kine¬ 
matically forbidden due to the low DM velocity. Tree- 
level spin-independent scattering no longer occurs, sig- 
nihcantly reducing all direct detection constraints. 

With spin-independent xi-nucleon cross section limits 
now in the vicinity of ^ 10“^® cm^ EllIIHl, other pro¬ 
cesses contributing to elastic DM scattering can also lead 
to constraints. These include higher-order loop-induced 
processes, as well as on-diagonal interactions suppressed 
by the mass splitting. Here, we give an overview of 
the various constraints, with details of the calculations 
provided in Appendix [B| 


Dark photon: Elastic xi-nucleus scattering can occur 
at loop level by exchanging pairs of dark photons via a 
virtual X 2 state (see Figs. 13 and 14). Whether elas¬ 
tic scattering enjoys a coherent nuclear enhancement de¬ 
pends on the characteristic energy scale of the dark pho¬ 
tons: the contribution peaks from loop momenta ~ 
and given a nuclear radius Rpf, a coherent enhancement 


-1 


occurs for m^' ^ Rjy suppressed for mx' ^ 

[ng. In this latter case, the DM resolves nuclear sub¬ 
structure and scatters predominantly off nucleons. In 
this paper, we consider mx' > > 100 MeV, and so 

we are nearly always in the regime where scattering oc¬ 


curs off of nucleons. Consequently, the elastic scattering 
rate is not only loop suppressed, but also suppressed by 
relative to conventional spin-independent scattering. 
For the parameter space of our study, we find that loop- 
induced direct detection is never important relative to 
other constraints on the dark photon model for splittings 
A/toi ^ 0.1, although direct detection could be more 
relevant for smaller splittings where collider constraints 
are less effective. 

In models where the Majorana masses for each DM 
charge eigenstate are maximally different — m^\ = 
mri + m^ in the notation of Sec. on-diagonal couplings 
to the dark photon are present as shown in Eq. (10). This 


gives rise to elastic scattering at leading order in e, but 
suppressed by the mass splitting A and the DM squared 
velocity. For mx' masses well below the Z, such that 
the dark photon mixes predominantly with the photon, 
the spin-independent cross section per-nucleon at zero 
momentum transfer in the limit of small splitting is 


lOTre^aao A^Z^ 


2 2 


(26) 


where v is the DM velocity and Z{A) is the atomic 
(mass) number of the target. We find that, due to 
the velocity suppression, the limits from LUX are 
subdominant to other constraints. For mA> ^ raz, 
an additional spin-dependent, velocity un-suppressed 
channel opens due to coupling with the SM axial current, 
but this constraint is also sub-dominant. 


MiDM: The dipole interaction induces elastic xi“ 
nucleus scattering at loop level (see Figs. and 15), 
and scattering can occur off of both nuclear charge and 
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Xi Xi 



FIG. 15: Leading elastic scattering diagram for direct de¬ 
tection in the MiDM scenario with the X 2 excited state in¬ 
tegrated out. The gray circle represents an insertion of the 
effective operator from Eq. (271. 


</ii <f>i 



FIG. 14: Leading elastic scattering diagram for scalar DM 
that couples off-diagonally to a vector mediator. Note that 
the A' vertex induces a third diagram at the same 

order. 


dipole [43l [74]. In the limit where the splitting is suf¬ 
ficiently large that X 2 cannot be produced on-shell, the 
leading effective interaction can be written in terms of a 
Rayleigh operator [30] 

( 4 ) . ( 27 ) 


where we have integrated out the excited state X 2 (we 
neglect additional velocity-suppressed terms). The direct 
detection scattering rate for this operator was computed 
in Ref. [30]. The resulting spin-independent, elastic cross 
section per nucleon is 


a 


(SI) 

n 


^2 0 J 


(28) 


where Z(A) is the atomic (mass) number of the target, 
Qo -- 76(0.3 -h 0.89y4i/3) fm'^ is the nuclear coherence 
scale, and is the DM-nucleon reduced mass. See 
Appendix [B] for a full discussion. We find that the LUX 
[64j constraint is subdominant over the parameter space. 


VI. CONCLUSIONS AND DISCUSSION 

In this article, we have proposed a set of searches at 
colliders that seek to exploit the distinctive signatures of 


inelastic dark matter (iDM). In iDM scenarios, direct and 
indirect detection constraints are weaker than for elastic 
DM, and therefore one prominent blind spot of current 
efforts is in the DM mass range of ~ 100 MeV—100 GeV. 
For fractional mass splittings between the two DM states 
of order ~ 10%, the heavier excited DM state can decay 
into the lighter state along with SM particles at decay 
lengths typically smaller than the size of collider exper¬ 
iments. These states can be reconstructed at colliders, 
giving rise to distinctive iDM signatures. 

These additional handles of iDM at colliders allow for 
targeted searches that may allow for sensitivity down to 
couplings favoured by a thermal DM origin in some of 
the simplest iDM scenarios. In particular, we focused on 
two representative models, described in Sec. |TTj that of a 
split DM system accompanied by a dark photon that ki- 
netically mixes with the SM, and that of magnetic dipole 
iDM. For each model, we propose new, powerful collider 
searches that can greatly increase sensitivity: in the dark 
photon model, we find that searches for monojet -I- miss¬ 
ing momentum -|- displaced leptons could cover nearly all 
of the remaining thermal relic territory for DM masses 
between ~ 100 MeV —100 GeV and splittings of 0(10%). 
For MiDM, the proposed searches for monojet -|- missing 
momentum -|- photon can also cover unexplored parame¬ 
ter space compatible with thermal freeze-out (for compa¬ 
rable masses and splittings as the dark photon model). 

In both scenarios, the thermal relic parameter space for 
much larger fractional splittings is essentially excluded 
by other experiments. For smaller splittings and DM 
masses below a few GeV, we find that monophoton -|- 
missing momentum searches at Belle II could provide 
complementary sensitivity in this regime |96j . Finally, 
we comment that significant strong sensitivity could be 
achieved for tudm 7 few GeV and small splittings with 
new electron beam-dump and active-target experiments 
[CTiiHTiiwi imfm] . proposed proton-beam fixed tar¬ 
get experiments |122Ill24j . and future direct detection 
experiments [3] optimized for low threshold sensitivity. 

Acknowledgments. We thank Ghris Brust, Ghristopher 
















20 


S. Hill, Bingxuan Liu, Maxim Pospelov, Philip Schuster, 
Abner Softer, Tim Tait, Natalia Toro, Wells Wulsin, and 
Itay Yavin for helpful conversations. The Perimeter In¬ 
stitute for Theoretical Physics is supported by the Gov¬ 
ernment of Canada through Industry Canada and by the 
Province of Ontario. El receives partial support from 
the Ministry of Research and Innovation - ERA (Early 
Research Awards) program. 


Appendix A: Freeze-out Rates 

In this Appendix, we summarize the relevant cross sec¬ 
tions and rates for the thermal relic abundances com¬ 
puted in Sec. |IV| We provide these rates to illustrate 
their characteristic scaling; consequently, we typically 
show only leading terms in an expansion in powers of 
DM velocity and mass splitting. For simplicity, we also 
present dark photon rates in the limit tuz, in 

which case we can assume that the dark photon couples 
predominantly to the electromagnetic current (this limit 
is also valid over much of the parameter space we con¬ 
sider). 


1. Annihilation Rates 


for the analytically intractable phase space for hadronic 
final states by using the known ratio of hadron and 
muon production in e“'“e“ annihilation R{s) = a{e^e~ —>■ 
hadrons)/cr(e+e“ —) from Ref. [T]. We decom¬ 
pose the full DM annihilation rate into SM states ac¬ 
cording to 


O'^’fuii = {av)e+e-S{s - Ami) 




9{s — Ami) + 0{s — Aml)R{s) 


• (A4) 


This is done in the vicinity of the hadronic resonances, 
^/s < 12 GeV. For larger DM masses, the couplings to 
partons are used. 

The R{s) approximation is valid only for couplings to 
the SM vector current; the SM axial current has a contri¬ 
bution to the cross section which differs from the vector 
current by terms proportional to where / is the fi¬ 
nal state fermion. In the dark photon model, the axial 
coupling only arises from mixing with the Z and scales 
like m\ilm\ for mA' ^ rnz] since the R[s) approxima¬ 
tion is only used for m^. < 6 GeV, and we consider mA' 
not parameterically larger than m^, the axial coupling is 
negligible and the R{s) approximation is valid. Similarly, 
for the dipole the Z contribution is suppressed by lm\ 
in the R{s) regime, and so is negligible. 


For the A' mediated interactions, the leading s-wave 
annihilation rate for fermion iDM is 


CT'c(XiX2-t //) 


Aire^aaoimi m.2)^ 
[(mi -I- m2)^ - m\,Y ’ 


where we have omitted corrections of order m //mi. For 
scalar iDM, the analogous A' mediated process is p-wave 
suppressed 


2. Scattering Rates 

For both fermion and scalar DM, scattering through a 
virtual A' mediator has the same parametric dependence 

T2 

(cru(x 2 e-t xie)) -- 167re^aaD^-. (A5) 

mA, 


(Tv{cl)i(j)2 -t //) 


8'!Te^aaumim2 
3[(mi -I- m 2 )^ — m^,]^ 


(A2) 


For fermion DM scattering through the dipole interac¬ 
tion, the leading order term in the velocity expansion is 


where v is the relative velocity. 

For the dipole mediated scenario, including both 7 and 
Z dipole operators, the s-channel annihilation rate into 
SM fermions // correct to first order in splitting and 
velocity is m 

ff) = aQ'jfil(l + 2vf — C[{mi+m2)'^] 

\ 

+ ivj + a})^('^[{mi+m2)‘^]) , (A3) 

/ 

where Q/ is the fermion electromagnetic charge, v{af) is 
the ratio of fermion’s vector (axial) charges to the elec¬ 
tromagnetic charge, and ({s) = s/(s — m|.). 

Over the DM mass window we consider in this pa¬ 
per (100 MeV— 100 GeV), annihilation to hadronic reso¬ 
nances can play a key role in setting the relic abundance 
in the early universe, yielding the spikes in the relic den¬ 
sity curves in Figs. [2][^ We can approximately account 


{av{x2e -t xie)) -- 3ap^. (A 6 ) 

In the above, the relative velocity u ^ 1 because the elec¬ 
trons remain relativistic. For typical values of the model 
parameters, these cross sections lead to rates larger than 
the Hubble expansion. 


3. Decay Rates 


For fermion DM in the A mediated scenario, the X 2 
width for A <?; mi, mA' is 


r(x2 -t X\ e+e ) 


157rm^, 


(A7) 


The corresponding expression for scalar DM de-exciting 
through a virtual A has the same value in the A <C mA' 
limit. 
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For excited dark-fermions decaying in the MiDM 
model, the width is 

r(x 2 ^ X 17 ) = —— ■ (A 8 ) 

TT 

In each scenario, the width is thermally averaged in 
the conventional way to obtain ( 72 )- 


parameterized by F{q) = (1 -|- and analogous 

form factors can be dehned as well for neutrons [12211126] . 
For the masses and splittings considered in this paper, 
existing direct detection constraints from LUX are too 
weak to constrain this process. 


2. A' Scalar iDM Direct Detection 


4. Inelastic Thermal Averaging 

To perform the thermal average for iDM in each sce¬ 
nario, we generalize the results of Ref. [125] to obtain to 
first order in the splitting A: 

W = dsa{s-so)y/sKi(^'^ , (A9) 

where sq = (wi -I- 77x2)^ and the averaging factor is 

N{T)=SmlmlTK^[^)K2['^) . (AlO) 

Ki and K 2 are the modified Bessel functions of the first 
and second kind. 


Appendix B: Direct Detection 


The scattering amplitude for this process is the sum 
of diagrams in Fig. 14 To leading order, the two box 
diagrams can be written in the HQET limit as [40| 


Ml -l- A.2 — 


g%e^t^{uu)n f d^\i\ 


Fi\e\y 


Anil J (27r)3 

d£0 (£ 0)2 {2mi-£/2f 

^ [(£ 0)2 - X {£)]3 


d£ 


(B3) 


where £ is a dimensionful Feynman parameter introduced 
to simplify denominator products of unequal mass di¬ 
mension, F{q) = {1 + q '^is the proton form factor 
introduced above, and A(£) = + m\, + \£ p. 

The third (triangle) diagram contributes 


AI 3 


2 2 2 /— \ 
gjje e (uu)n 


r.c [ rde° F(|£|)2 

Jo y ( 2^)3 1^277 [(€ 0)2 _a(£)] 3 ’ 


1. A' Fermion iDM Direct Detection 


SO the total amplitude A = i® squared to give 


At loop level, iDM mediated by a massive A can scat¬ 
ter at dire ct detection experiments via the box diagram 
in Fig. |l3[ For ntA' 'A the 1-loop Feynman diagram 
in Fig. |13| involves a single nucleon on the SM side (as 
opposed to a nucleus). The differential cross section can 
be defined as 


da 

d£l 


Att^ 


[/(m^, A)| 


(Bl) 


where / is the integral over the loop momentum and 
is the y-nucleon reduced mass. This can be evaluated 
using the formalism of Heavy Quark Effective Theory 
(HQET) as in Ref. [10] . For zero momentum transfer, 
the integral can be written as 


/(tox. A) 


Sie'^aaD 

A 



d\m^ 


Fm? 

(|€|2 + m^,) 2 ’ 


(B2) 


for simplicity, we show the result only for the A cou¬ 
pling to the electromagnetic current, which is valid in 
the 777(4 ^ limit. 

In contrast with Refs. gomini, the A in our regime 
of interest generates a short range force that resolves nu¬ 
clear substructure, so the scattering is predominantly 
off nucleons without the additional enhancement 
from coherent, loop-scattering off nuclear targets. The 
monopole form-factor for scattering off protons can be 


i\A\^)=Aml\I(jni,A)\\ 
where the integral is 


(B4) 


/(ttii. A) = - 


27771 


Siaajje^ 

Ami 


d\e\\£f 


Am\ 


+ 


(|€|2+TO^,)3/2 |£| 2 +^ 2 ^y 4(|€|2+m^,)2 


|€|2 + m^,)" 

Attii I 


Fmf 

(B5) 


and we have averaged over the spin of the nucleon target. 
The 01 -nucleon cross section is now 


_ |/(777i, A)|2m^ 
477(7771 + 777„)2 


(B 6 ) 


3. Magnetic Dipole Direct Detection 

For magnetic dipole scattering through the same boxes 
in Fig. |15[ we can work in the limit where A is much 
larger than the kinetic energy of the DM, and we in¬ 
tegrate out the excited state. This is equivalent to per¬ 
forming an HQET expansion on the DM propagator. The 
leading-order effective interactions between DM and nu¬ 
clei arises from two photon interactions, which at leading 
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order in the expansion and up to Lorentz index contrac¬ 
tions are 

/:= + tXi7^XiF^.F^^) , (B7) 


where qis the momentum transfer, = mimpf /{mi + 

mpp) is the DM-nucleus reduced mass, and Qo = 
76(0.3is the coherence scale for a target 
of mass number A. In this limit the spin-independent 
cross-section per nucleon can be approximated as 


where = (l/2)e^i,„^F“^ is the dual electromagnetic 
field-strength tensor. To match the conventions in [ID] for 
this limit, we identify = /r^/2m2. At one loop 

level, these interaction give rise to elastic scattering at 
one loop level through the diagram depicted in Fig. [Ts] 
however, the contribution through scattering from the 
second operator in Eq. (B7) is velocity suppressed and 


can be neglected when computing direct detection rates. 

The differential cross section for coherent non- 
relativistic scattering off nuclear target N is 


da 

dcosO 


l^xN 

2tt 


M7Q0 f\q\^ 

2m2 ^ V Qo 


(B8) 


a 


(SI) 

n 




(B9) 


where /i^„ DM-nucleon reduced mass and and we have 
approximated the Helm form factor |127] as Fn{0) ~ 
2/771 at zero momentum transfer. Note that this form is 
only valid in the limit where the mass splitting is much 
larger than the DM kinetic energy, which is trivially sat¬ 
isfied in our regime of interest throughout this paper. 
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